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Abstract. In this paper we explore the possibility of using computa- 
tional algebraic methods to analyze a class of consensus protocols. We 
state some necessary conditions for convergence under consensus proto- 
cols that are polynomials. 

Key words: Grobner Bases, Dependency graph, Algebraic Variety 

u 

**i ' 1 Introduction 

CZ2 

Consensus among agents is an important problem in multi-agent systems and 
congestion/flow control in communication networks. We formulate the problem 
as follows. 

Consensus protocol is an ordered set of functions (/i,/2, • • • , /jv)j where 
qq ' fi : R N — > R, i = 1, ...,JV. We consider the following first order consensus 

iy) , dynamics among n agents under this protocol 

m 
c^ 

Here xi,X2, ■ ■ ■ , xn are the states of agents, the initial values of which are 
assumed to be known to the respective agents. The problem of interest is to 
determine whether the system of equations in ([1]) attains a consensus. Consensus 
is an equilibrium at which x\ = X% = . . . = x^ . In other words, we are interested 
K^ . to know whether there exist an a £ 1 such that when x = x e , where x = 

$— i ' (x±, X2, ■ ■ ■ ,%n), and x e is of the form x e = (1, 1, . . . , l)a, we have /i(x e ) = 

0,i — 1, . . . , N, and starting from an arbitrary x, the system ([T]) will converge 
to an equilibrium solution x e . Equation [T] can be written compactly as 

x = f (x). (2) 

where f = (/i, / 2 , . . . , /at)- Here f , the ordered AT-tuple, is the consensus protocol 
and we are interested in protocols which will take the dynamical system ([1} to 
consensus. 



fi(Xl,X2,-.-,XN) i = l,...,N. (1) 
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It is possible to associate a graph with every dynamical system of the form 
(fTJ). Let Q = (V, £) be a such graph with vertices V and edges £. Then V = 
{1, 2, ... , N} and a directed edge e^ = (i, j) from node i to node j belongs to 
£ if /j is a function of a^. If fj is a function of Xi, we say that Xi belongs to the 
support of fj, that is, xt € supp(fj). Following Stigler [T], we call the graph Q 
the dependency graph of (fl]). The neighbors of an agent /node x% are all those 
nodes which are in the support of /,. Thus, the set of all neighbors of a node i is 

Afi = {j : eji £ £}. (3) 

The most widely analyzed protocol for consensus [5] is a linear protocol of form 

ii = ^2 (xj - x^, i = l,,..,N. (4) 

The adjacency matrix, A — [o»j], associated with a graph Q = (V, £) is defined 



as 

1 if eij G £ 

otherwise. 

A diagonal matrix called the degree matrix, D = [da], is defined as da — 
number of edges directed towards Xi, which is equal to Ylj=i a ij- r ^^ ie g ra P n 
Laplacian, L, is defined as L = D — A. Using the definition of graph Laplacian, 
we can rewrite ([4]) as 

x = -Lx (5) 

There is a great amount of literature on consensus under linear protocols. The 
assignment of a graph structure to the problem makes the problem amenable 
to application of tools from algebraic graph theory like graph Laplacian [3], 
stochasticity of nonnegative adjacency matrices |3], etc. to analyze the consensus 
in system given by ^ or a normalized form of it. However, these methods for 
analyzing stability and convergence to consensus of the system of equations may 
not be useful for nonlinear protocols. 

There has been a few papers, in consensus literature, addressing consensus 
under nonlinear protocols. Olfati-Saber and Murray [3J considers consensus un- 
der a nonlinear protocol of form 

Xj = y^ hj(xj - Xj), i = l,...,n (6) 

where hi,i = l,...,n are functions which are uneven, locally Lipschitz, and 
strictly increasing. 

Liu and Chen [5] extends this to protocols of form 

Xi = J2 ( h (Xj) ~ Hxi)) , i = l,-..,N (7) 

3&/V 4 

where h is assumed to be an increasing function. However, this protocol demands 
that h be same for all agents. 
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Moreau [H] proves the convergence of consensus protocol under the assump- 
tion that the protocol is such that, in the discretized version of it, the updated 
value of a node is a strict convex combination of current values of the node and 
its neighbors . This implies that the continuous time protocol of the form 

Xj = y^ hj(-*)(xj —Xj), i = l,...,N (8) 

will achieve consensus. 

This paper we pose problem of consensus under polynomial protocols in the 
framework of computational algebra and give necessary conditions for the conver- 
gence. This paper is organized as follows. In §[5] we give preliminary observations 
on polynomial consensus and review basic background in computational algebra 
that is required. We present main results in § [3] and give concluding remarks in 
§11 

2 Preliminary Observations and Basic Computational 
Algebra 

2.1 Preliminary Observations 

We give few necessary conditions for convergence to a consensus under polyno- 
mial protocols, using the tools and language of algebraic geometry. It is clear 
that for a consensus to be achieved, the protocol f should use all the nodal val- 
ues. A consensus is not achieved, except for some particular initial conditions, if 
values of one or mode nodes are not used in the dynamics © . This leads us to 
the following proposition. 

Proposition 1. If the system ii — fi,i = 1,...,N, achieves consensus, then 
for all j 6 {i, . . . , N}, there exists i S {i, . . . , N} such that Xj G supp(fi). 

The proposition asserts that the dependency graph of a protocol that achieve 
consensus should not have any isolated nodes. 

Also, if more than one nodal values are not being updated, it is not possible 
to arrive at a consensus in general. However, there might exist a few particular 
initial conditions for which a consensus on nodal values is reached. For example, 
the initial values of the nodes whose values are not being updated are equal 
and equal to the consensus value achieved by other nodes. However, this is 
an exception rather than a rule. Therefore, we have the following necessary 
condition. 

Proposition 2. If more than one ft are zero polynomials, then the system ii = 
fi,i = 1, . . . , N , will not achieve consensus. 

A result from theory of linear consensus protocols states that a consensus is 
not achieved if the underlying graph does not have a directed spanning tree [7] . 
Proposition [5] has the same flavor. 
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From here onwards when we say that a protocol f satisfies the conditions for 
Propositions [1] and [U we mean that the protocol involves all the nodal values 
and at most one of the polynomials in f is a zero polynomial. And, whenever 
one of the polynomials is a zero polynomial, by the set {/i, . . . , /jv} we mean 
the set without the zero polynomial. 

If the system of equations in Eq. ([2]) leads to a consensus, then we have x = 0, 
which implies f(x) = 0, asymptotically. The equilibrium points of the system of 
equations in Eq. ^ are the 'roots' of the set of polynomial equations f (x) = 0. 
Now we need some notions of algebraic geometry and Grobner bases. 

2.2 Basics of Algebraic Geometry and Grobner bases 

Throughout this paper, k represents a field (e.g., R, C). Set of all monomials 
in indeterminates x\, . . . ,xm is denoted by Z™ and set of all polynomials in 
indeterminates x\, . . . , xm with coefficients in k is denoted by k[x\, . . . , xjq]. Let 
/i, . . . , f s e k[xi, . . . , xjv]- We use the notation V(/i, . . . , f s ) — V to represent 
the varieties, where 

V = {(ci, ...c n )ek N : fi( Cl , ...c N ) = 0,1 <i<s} . 

V is uniquely determined by the ideal generated by /i , . . . , f s . This ideal is 
denoted by (f%, . . . , f a ) and hence we have 

V(/i,...,/ s ) = V«/i,...,/ s » . 

The set of equilibrium points of most of the consensus protocols contain the 
subspace x\ = ■ ■ ■ = x^. Let us call this subspace where states of all agents 
are equal as S. If a protocol f leads to consensus, then it should have at least 
one equilibrium point belonging to the subspace S. Thus we get a necessary 
condition for convergence of Eq. ([2]) as follows. If Eq. ([2]) converges, then SO V(f ) 
is nonempty. Geometrically, the subspace S C k N is a solution to set of equations 

s = {xi-x 2 ,x 2 -x 3 ,...,x N -i-x N } (9) 

In the language of algebraic geometry, S is the variety of s, that is, S = V(s). 
Therefore, a consensus would imply V(s) n V(f ) is non-empty. 

Let J be the ideal generated by set of polynomials s describing the subspace 

5 as described before. It is easy to see that the affine variety of the set of 
polynomials s is equal to the affine variety of the ideal J generated by them [5] , 
that is, V(J) = V(s) = S. 

Since S is an irreducible variety (cannot be written as the union of non-empty 
varieties), the ideal J is prime [8:. Since J is prime, the radical ideao of J is 
itself. Now, by Hilbert's Strong Nullstellansatz [S], the ideal 2T(5), the ideal of 
all polynomials that vanishes at every point on S, is J. 



The radical ideal of J, denoted by \fj, is the set {/ £ 1l\xi,...,Xn] ■ f" £ 
J for some n G Z>o} 
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Throughout this article, we will be using some results from algebraic geom- 
etry. These will be given as propositions and theorems without proofs and we 
refer an interested reader for details to Cox et al. [8]. 

Theorem 1. f^ If I and J are ideals in k[xi, . . . ,xn], then V(I + J) = V(I) fl 
V(J). 

Here, I + J is defined as follows. 

Definition 1. If I and J are ideals of ring k[xi, . . . ,xn], then the sum of I 
and J, denoted by I + J , is the set 

I + J = {f + g\feIandgeJ}. 

In fact, / + J is the smallest ideal containing both / and J [8]. We have the 
following proposition. 

Proposition 3. Let I be the ideal generated by {/i, . . . , /jv}, fi G k[xi, • ■ • , Xn], 
i = 1, . . . , N and J be the ideal (xi — X2, x-i — X3, . . . , xpj-i — Xpf). Let S C k 
be the subspace given as S = {(xi, . . . , xjv) G k w : Xi = X2 = • • • = x^}. Then, 
if the polynomial dynamical system &i — fi, i — 1, . . . , N attains consensus, then 
V(I + J) is not empty and I + J is a proper subset of k[xi, . . . , xn] ■ 

Proof. A necessary condition for consensus is that at least one of the stationary 
points of the system Xj = /,, i = 1, . . . , n should belong to iS. Set of stationary 
points or the invariant set of this system is the solution to set of polynomial 
equations f\ = 0, . . . , /at = which is given by V(I). Thus the necessary condi- 
tion demands that V(I) PI S is not empty. Since S is Zariski closeqj as S is the 
variety of ideal J , we get V(I) fl V(J) ^ 0. This along with Theorem Q] gives 
that V(I + J) is not empty. Since V(I + J) is not empty, by Hilbert's Weak 
Nullstellensatz [S], / + J is proper ideal of k[xi, . . . , xjv]. 

The statements on ideals can be checked by calculating their Grdbner bases 
which can be done using any symbolic package that supports algebraic geometric 
calculations, given the set of polynomials that generate the ideal. Now, we define 
a Grobner basis. 

Grobner bases is a generalization of the division algorithm in a single variable 
case (k[x]) to the multivariate case (k[xi, . . . , Xjv]). Grobner basis also generalizes 
Gaussian elimination in linear polynomials to nonlinear polynomials. For the 
multivariate division algorithm, we need the notion of monomial order. 

Definition 2. A monomial order or term order on k[x\, . . . ,xn] is a re- 
lation ~< on Z> that satisfies following conditions (i) -< is a total ordering on 
Z> , (ii) if a < fi, for a, /3 s Z> then for any 7 £ Z^ it holds a + 7 -< f3 + j, 
and (Hi) -< is a well-ordering on Z> . 



5 A set belonging to k^ is said to be Zariski closed if it is the solution to a set of 
polynomials in k[a;i, . . . , xn]- 
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Given such an ordering -<, one can define the leading term of non-zero polynomial 
/ £ k[xi, . . . , xn] as a term of / (the coefficient times its monomial) whose 
monomial is maximal for -<. We denote this leading term by LT^(/) and the 
corresponding monomial by LM_/(/)- 

Definition 3. An ideal a C k[xi, . . . , xn] is said to be a monomial ideal if there 
is a set A C Z> , possibly infinite, such that a = (x a : a <E A). 

Given any ideal a C k[xi, . . . ,Xn], the ideal defined as (LM^(/) : / £ a) is a 
monomial ideal and is denoted by LM^(o), which is known as leading monomial 
ideal of o. By Dickson's lemma [8j p. 69] , the ideal LM^ (a) is generated by a finite 
set of monomials. Dickson's lemma and the multivariate division algorithm leads 
to a proof of Hilbert bases theorem which states that every polynomial ideal can 
be finitely generated, which further lead to a definition of Grobner basis [5J § 2.5]. 

Definition 4. Fix a monomial order -< on k[xi, . . . , Xjv]- A finite subset G = 
{gi, . . . , g s } of an ideal a C k[xi, . . . , xjv] is a Grobner basis if and only if 

LM,(o)'=(IM,( 5l ),...,I¥,y). 

Given a set of generators of an ideal, the Buchberger algorithm [TU] can be used 
to compute a Grobner basis of the ideal with respect to various term orders. 
The algorithm and its variants are implemented in most symbolic computation 
programs. Note that a Grobner basis is not unique, but one can transform it to 
a reduced Grobner basis which is unique for every ideal in k[xi, . . . , xjv]- In the 
sequel, when we say 'the Grobner basis', we mean the reduced Grobner basis. 
The Buchberger algorithm provides a common generalization of the Euclidean 
division algorithm and the Gaussian elimination algorithm to multivariate poly- 
nomial rings. 

3 Consensus under polynomial protocols: Necessary 
conditions 

Now, an equivalent statement of Proposition [3] in terms of Grobner basis can be 
given as follows. 

Corollary 1. If G is a Grobner basis of I + J where I and J are as defined in 
Proposition [3J then the affine variety of G is not empty and the ideal generated 
by G is a proper ideal of k[xi , . . . , xn] ■ 

The above condition will always be satisfied if none of the fi,i — 1, . . . , N, has 
a constant term, as then, x = is a common solution of polynomials in both / 
and J. We call this the trivial consensus where a consensus occurs if the initial 
conditions on all the nodes are simultaneously zeros. Similarly, consensus occurs 
for other special initial cases too. Since the conditions of Proposition[3jalso holds 
for protocols that can achieve trivial consensus also, it is a weak result. 

From Proposition [3] we have that, to show that a system does not attain con- 
sensus, it is enough to show that I+J is the entire polynomial ring k[xi, . . . , xjv], 
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or equivalently by Hilbert's Weak Nullstellansatz, the Grobner basis of ideal I+J 
is 1 (since 1 generates the whole ring k[xi, . . . , Xjv] which has an empty variety). 
For this, we need a way to calculate the ideal I + J given / and J, which the 
following proposition gives. 

Proposition 4. ' L 8llf I = (/i, • ■ • ,/r) an d J — (gi, ■ ■ ■ ,g s ), then I + J = 
(fi,---,fr,gi,---,9s)- 

Thus to check the Proposition^ given (f), the ideal generated by protocols, and 
(s), the ideal the variety of which is S, we need to calculate the variety of ideal 
(f ) + (s) which by Proposition @] is equal to the variety of (f , s). Since a Grobner 
basis of (f, s) has the same variety as that of (f, s), it is enough to calculate 
the Grobner basis of (f , s) and look at its variety. This can be done using any 
software that support computation of Grobner basis. For the examples given in 
this chapter, we use Mathematica. 

Example 1. We consider a nonlinear protocol as follows. 

±i = x\ + x 2 + 1, ±2 — Xl + x 3 + 3, ± 3 = x 2 - x 3 . 

GroebnerBasis[{{xi 4- X2 + 1, x\ + X3 + 3, X2 — x$},{x\ — X2, X2 — X3}}] 

{1} 

For the system in Eq. ([1} , a consensus is achieved if there exists an a £ R 
such that x\ = x 2 = ■ ■ ■ = Xn = a is an equilibrium point. However, systems 
for which X\ = x% = • • • = xjq = a is an equilibrium point for all a £ R 
are of particular interest in the consensus community. Therefore, in the rest of 
the paper, we consider only such protocols. The necessary condition as given 
in Proposition [3] then amounts to V(J) C V(I), where I = (fi, / 2 , . . . , Jn) and 
J = (x 2 -xx,x 3 -Xh,...,Xn-i -%n)- 

A consensus cannot be achieved over a graph that is disconnected. Thus a 
necessary condition for a consensus to occur under a protocol is that its de- 
pendency graph should be connected. In the case of a linear protocol with a 
dependency graph that is directed, it should be strongly connected (or at least 
should have a directed spanning tree) for consensus to occur [7]. A result from 
algebraic graph theory states that if a directed graph is strongly connected, then 
for a linear consensus protocol the graph Laplacian matrix associated with it is 
irreducible |11) . This essentially means that by a permutation or an automor- 
phism of the variables, the associated graph Laplacian can be made to be of 
a block diagonal form. The number of blocks are equal to the number of dis- 
connected components of the graph [7J. Given a protocol f = {/1, ■ • ■ , /jv}, its 
dependency graph Q = (V, £) can be found. It is possible to characterize the 
properties of this graph by defining the following relation. 

Definition 5. Given two nodes i,j £ V, we say i ^ p j if there exists a path from 
node j to node i. If i -< p j and j -< p i, then nodes i and j are path- equivalent 
and we denote this by i ~ p j. 
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The relation -< p is anti-symmetric since if i -< p j and j -< p i, then i ^ p j. 
It is transitive. The above relation is reflexive and hence a partial order if the 
dependency graph contains self-loops. 

We have the following results that are immediate from the definition of the 
order relation -< p . 

Proposition 5. Let Q = (V, £ ) be graph with V — {1, . . . , N}. Let Maximal^ V 
denote the set of maximal elements in V under the order relation -< p . Then, 

1. The graph Q is strongly connected if and only if Maximal^ V = V 

2. The graph Q has a directed spanning tree if and only if 
# {Maximally j ~ p ) = 1 

where # denotes the cardinality of a set. 

The necessary and sufficient condition for achieving a consensus under linear 
protocol over a graph is the existence of a directed spanning tree [?]• In fact, the 
existence of a directed spanning tree is a necessary condition not only for the 
linear consensus protocol but also for any protocol that achieves consensus over 
a graph. 




Fig. 1. Dependency graph with a directed spanning tree rooted at node 3. 




Fig. 2. Consensus under protocol in Eq. (|TU|) . 
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Example 2. The linear protocol 






IMI 



applied to dependency graph given in Fig. [T]is 



x\ = x 2 
X4 = x 2 



x 5 - 

X4,, 



2*1, ±2 = Xl 
X 5 = Xl 



%3 



X 4 — 3x 2 , X3 = 0, 
2x 5 . 



Convergence to consensus under this protocol is given in Fig. [21 
Example 3. We now consider consensus under a nonlinear protocol 



E 



I M I 



jeM 



for the dependency graph in Fig. [T] 



X4 Xn 

_ 2 
X 4 — X o 



X5 — 2x1, X2 = x 



Jb ^ <*J I I - tL Q 



X4 ^^2 ) *^3 

2x^. 



0, 



Figure [3] shows the convergence to consensus of this protocol. 



(10) 



(11) 



Fig. 3. Consensus under protocol in Eq. (|TT|) . 



Example 4- The third protocol considered earlier is 



n 



„IM| 



Applied to dependency graph of Fig. \T\ we get 



(12) 



Xl = X2X5 — x\, X2 = X1X3X4 — x\, ±3 =0, X4 = X2 — X4, X5 = X1X3 — x\. 



Figure H] shows that this protocol leads to consensus. 
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Fig. 4. Consensus under protocol in Eq. (|T^|) . 



Examples HH1] give protocols all of which had a dependency graph that is not 
strongly connected but still achieve consensus. This shows that the dependency 
graph of a protocol being strongly connected is not a necessary condition for 
consensus to occur. 

Using the elimination theorem from algebraic geometry, we can get a result 
that enables us to determine whether the dependency graph of a polynomial 
protocol has a directed spanning tree. To state the elimination theorem, we 
need concepts of elimination order and elimination ideal as defined below. 

Definition 6. Consider k[xi, . . . , x n , yi, . . . , y m ] a polynomial ring in indeter- 
minates x\, . . . , x n , y\, . . . , y m . Let -< x and -< y be monomial orderings on x and 
y variables respectively. Define an ordering relation ~^[{x}y{y}] on ^>^ m (i-e set 
of all monomials in indeterminates Xi, . . . , X n , t/i, . . . , y m ) as follows: 



x a y b -<[{*»-{»}] xC v d 




c and b -< y d 



where a,c £ Z> and b,d € Z> . The term order "^[{x}y{y}] * s called elimina- 
tion order with the x variables larger than the y variables (which is indeed a 
term order). 

Definition 7. Given I = (/i,...,/ 8 ) C k[xi, . . . , xjv] the l-th elimination 

ideal // is the ideal of k.[xi + i , . . . , xn] defined by Ii = I n k[xj+i , . . . , xn] 

Now, we sate the elimination theorem. 

Theorem 2. ' L 8l Let L C k[xi, . . . , xm] be an ideal and let G be a Grobner Basis 
of L with respect to an elimination order where {x\, X%, . . . , x{\ >- {iC|+i, • • • , xn} 
for every < I < N . Then the set 



Gi = Gnk[j|+i, . . . ,x N ] 
is a Grobner Basis of the l-th elimination ideal I{. 
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If the dependency graph of a linear consensus protocol contains a directed span- 
ning tree, then the corresponding Laplacian matrix, L, will be of rank JV-1 0- 
This means that it is possible to eliminate N — 2 variables from the protocol 
Lx using Gaussian elimination. Since the elimination in the multivariate poly- 
nomial case is the generalization of the Gaussian elimination, we expect this nice 
property of being able to eliminate N — 2 variables in the linear protocol case 
to carry over to polynomial protocols that attain consensus. 

Theorem 3. Let ft E k[xi, . . . , x N ], i = 1, ...,JV. Let I = (/i,/ 2 , . . . ,/jv) 
be such that V(J) C V(I) for J = (x2 — Xi,X3 — x, . . . ,Xjv-i — %n)- If the 
system Xj — ft,i = 1,...,N attains consensus, then the (N — 2)-th elim- 
ination ideal of (/i, • • • , /at) is not empty under any elimination order with 
{x a m, . . . ,aV(Ar-2)} ~> { x a{N-i), x (j(N)} where a belongs to the set of every 
possible permutation of {1, . . . , N}. 

This follows from the fact that a consensus is achieved only if the dependency 
graph has a directed spanning tree. We illustrate this through the following 
example. 

Example 5. Consider a polynomial protocol given as 

Xl=X2—Xl, X2=Xl—X2, ±3 = 0, ±4 = XiX 2 X 3 - X% 

The corresponding dependency graph is shown in Fig. [5] This dependency graph 
does not have a directed spanning tree and thus will not attain a consensus. 
Thus by Theorem [31 there should exist an elimination order under which it may 
not possible to eliminate 2 variables from the Grobner Basis of {x 2 — X\,x\ — 
X2, £i£2:£3 — x\}. If we try to eliminate the highest two variables under the 
elimination order {x^x%} -< {x 2 ,xi}, we get an empty set implying that the 
elimination is not possible. 

GroebnerBasis[{x2 — xi,xi — X2,xiX2X3 — x|}, {x4,X3}, {x2,xi}] 
{} 




Fig. 5. Dependency graph of protocol in Example [S] 



It is possible to extend the result in Theorem [31 to the following. 

Proposition 6. Let f = (/i,--.,/jv) S k[xi,. . .,Xn] N be such that V(J) C 
V(I) for J = (x 2 - Xi,x 3 - x2, ...,xjv-i - Xjv) and L = (/i,/ 2 , . . . ,/jv)- Let 
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I G N be the smallest I such that the l-th elimination ideal of (/i,-.-,/jv) * s 
empty under an elimination order with {xi x , ■ ■ ■ , Xi,} >- {Xi l+1 , . . . ,Xi N } where 
{ii, . . . ,%n} is a permutation of {1, ... , N}. Then, # yMaximal^ p V / ~ p ) = N — 
I. 

For Example [5l we have N = 4, I = 2, and from Fig. [5] we see that 

# (Maximally/ ~ p ) = 2. 

4 Concluding Remarks 

We explored a novel way of looking at the consensus over networks-the algebraic 
geometric way-when the protocols are polynomials. We gave several necessary 
conditions for obtaining consensus under polynomial protocols. This enables one 
to comment on the convergence to consensus of a protocol by looking at the 
computed Grobner Basis of the set of polynomials in the protocol. We also gave 
a sufficient condition for convergence under polynomial consensus. We conclude 
by remarking that algebraic geometry, if used properly, has sufficient powerful 
tools to analyze consensus under switching protocols and to design consensus 
protocols with desired behaviors. 
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